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Abstract We present a data-adaptive spectral method - Monte Carlo Singular Spec¬ 
trum Analysis (MC-SSA) - and its modification to tackle astrophysical problems. 
Through numerical simulations we show the ability of the MC-SSA in dealing with 
1//^ power-law noise affected by photon counting statistics. Such noise process is 
simulated by a first-order autoregressive, AR(1) process corrupted by intrinsic Pois¬ 
son noise. In doing so, we statistically estimate a basic stochastic variation of the 
source and the corresponding fluctuations due to the quantum nature of light. In ad¬ 
dition, MC-SSA test retains its effectiveness even when a significant percentage of 
the signal falls below a certain level of detection, e.g., caused by the instrument sen¬ 
sitivity. The parsimonious approach presented here may be broadly applied, from 
the search for extrasolar planets to the extraction of low-intensity coherent phenom¬ 
ena probably hidden in high energy transients. 
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Singulai' Spectrum Analysis for astronomical time series 


1 Colored noise and MC-SSA 

1//^ power-law noise is known to be highly relevant in several astrophysical sys¬ 
tems [4, 5, 3]. It includes the well-known white noise (j3 = 0), pink noise (j3 = 1) 
and red or Brownian noise (j3 = 2). Such kind of noises are very often not tied to 
instrumental disturbance but rather they are property of the observed sources that 
emit radiation varying in a stochastic manner. In the context of our work, it is suffi¬ 
cient to qualitatively characterize 1 //^ to generate a parsimonious hypothesis test in 
MC-SSA (see below). Subsequently, the dynamical origin of noise fluctuation or the 
extra-noise variance - according to acceptance or rejection of our null hypothesis - 
has to be determined by theoretical analysis. 

Singular Spectrum Analysis (SS A) is an effective, data-adaptive and non-parametric 
method for the decomposition of a time series into a well-defined set of indepen¬ 
dent and interpretable components that include a non-linear trend, anharmonic, 
amplitude-modulated oscillations, and noise [6]. In its later developments, the 
Monte Carlo approach to signal-to-noise separation introduced by [1] has become 
known as Monte Carlo SSA (MC-SSA) [2]. The distinction between what is noise 
and what is signal is made through the measure of the resemblance of a given noise 
surrogate to the original data via eigendecomposition of the time-lagged covariance 
matrix. For this study the noise surrogates set is generated by using an AR(1) pro¬ 
cess. In practice, the AR(1) coefficients are estimated from the time-series under 
consideration by using a maximum-likelihood criterion. Subsequently, these AR(1) 
noise surrogates are corrupted by Poisson noise to mimic the effect of a photon 
counting detector. 


2 Simulated 1 //^ series test 


The procedure has been tested on a large sample of artificial time series of arbitrary 
colored noises. As a test of validity, we show several discrete colored noise vectors 
of length N = 5000, with a power-law distribution of slope p ranging from 1 to 
2 in steps of two. To obtain sequences of colored noises we use the Matlab library 
CNOISE' . A threshold of signal identification is incorporated by subtracting a mean 
value ~ 35% of the maximum peak-flux, C,nax- In this way, we take into account 
episodes of emission during which the emission count rate drops to the background 
level. Finally, the time series are corrupted by Poisson noise to simulate the effect 
of the shot noise based on Poisson photon statistics. The resulting test series, are 
shown in the supplementary video^. The signal to noise ratio (SNR) of the artificial 
series test increases as the simulations run. 


* http://people.sc.fsu.edu/ jburkardt/m_src/cnoise/cnoise.html 
^ https://vimeo.eom/120699083 
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3 Conclusion and future work 

We present a modified data-adaptive spectral method - Monte Carlo Singular Spec¬ 
trum Analysis (MC-SSA) - in which a parsimonious (dynamic and instrumental) 
noise-model Hq is adopted for astrophysical applications, i.e. Hq = {AR(1) + Pois¬ 
son noise}. AR(1) noise takes into account the long-term variability in the power- 
law slope -p of a light-curve’s spectral density, and the Poisson noise considers the 
short-term variability imposed by the quantum nature of light. Our simulations show 
a remarkable effectiveness of the model for colored noises with j3 values between 
1.5 and 2 and with Qnax > 1000 counts; otherwise one should re-bin the light-curve 
to improve the signal-to-noise ratio of the source. 

All analyses were performed by using SSA-MTM Toolkit freeware that has been 
developed by the Theoretical Climate Dynamics Group at UCLA^. To obtain pre¬ 
sented results we relied on an advanced option in SSA-MTM Toolkit for MC-SSA 
that allows to read in ensemble of surrogates created off-line and corresponding to 
an arbitrary noise null-hypothesis Hq. In this regard, MC-SSA allows to iteratively 
refine the working hypothesis and adapt the null-hypothesis Hq for various theoret¬ 
ical scenarios. We plan to extend further SSA-MTM Toolkit for astrophysical ap¬ 
plications by directly including the algorithm for parsimonious Hq in the freeware. 
Current and future work will focus on SSA as a promising technique for analysing 
complex spatio-temporal structures to detect possible periodicities and other statis¬ 
tical regularities due to various astrophysical processes. 
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